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ABSTRACT 

Context. The evolved main-sequence Sun-like stars KIC 10273246 (F-type) and KIC 10920273 (G-type) were observed with the 
NASA Kepler satellite for approximately ten months with a duty cycle in excess of 90%. Such continuous and long observations are 
unprecedented for solar-type stars other than the Sun. 

Aims. We aimed mainly at extracting estimates of p-mode frequencies - as well as of other individual mode parameters - from the 
power spectra of the light curves of both stars, thus providing scope for a full seismic characterization. 

Methods. The light curves were corrected for instrumental effects in a manner independent of the Kepler Science Pipeline. Estimation 
of individual mode parameters was based both on the maximization of the likelihood of a model describing the power spectrum and 
on a classic prewhitening method. Finally, we employed a procedure for selecting frequency lists to be used in stellar modeling. 
Results. A total of 30 and 21 modes of degree 1=0, 1,2- spanning at least eight radial orders - have been identified for KIC 10273246 
and KIC 10920273, respectively. Two avoided crossings (/ = 1 ridge) have been identified for KIC 10273246, whereas one avoided 
crossing plus another likely one have been identified for KIC 10920273. Good agreement is found between observed and predicted 
mode amplitudes for the F-type star KIC 10273246, based on a revised scaling relation. Estimates are given of the rotational periods, 
the parameters describing stellar granulation and the global asteroseismic parameters Av and v max . 

Key words, methods: data analysis - stars: individual (KIC 10273246, KIC 10920273) - stars: solar-type - stars: oscillations 



1. Introduction 

The NASA Kepler Mission was designed to use the transit 
method to detect Earth-like planets in an d near the ha bitable 
zones of late-type main-sequence stars dBorucki et alj l2010t 



* e-mail: campante@ast.ro . up . pt , campante@phys .au.dk 



iKochet al.ll20Toh . The satellite consists of a 0.95-meter aper- 
ture telescope with a CCD array and is capable of producing 
photometric observations with a precision of a few parts-per- 
million (ppm) during a period of 4-6 years. The high-quality 
data provided by Kepler are also well suited for conducting as- 
teroseismic studies of stars as part of the Kepler Asteroseismic 
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Investigation (KAI; iGilliland etaD 1201 Oal) . Photometry of the 
vast majority of these stars is conducted at long cadence (29.4 
minutes), while a revolving selection of up to 5 12 stars are moni- 
tored at short cadence (58.85 seconds). Short-cadence data allow 
us to investigate solar-like oscillations in main-sequence stars 
and subgiants, whose dominant periods are of the order of sev- 
eral minutes dChaplin et al]l2010t IChristensen-D alsgaard et all 
l2010tlMetcalfe et al.ll20l6Th 

The information contained in solar-like oscillations 
allows fundamental stellar properties (e.g., mass, ra- 
dius and age) to be determined (e.g.. | Stello et at] 12 009b; 
IChristensen-Dalsgaard et alJl2010t iKallinger et al.ll2010ah . The 
internal stellar structure can be constrained to unprecedented 
levels, provide d that individual o scillation mode parameters are 
extracted (e.g., Cunha et al. 2007). This is possible in the case of 
the highest signal-to-noise ratio (SNR) observations, leading us 
to hope that asteroseismology will produce significant improve- 
ment on the theories related to stellar structure and evolution, 
on topics as diverse as energy generation and transport, ro tation 
and stellar cycles (e.gjKaroff et al.l l2009fc iGarcia et al . 2010). 

Solar-like oscillations in a few tens of main-sequence 
stars and subgiants have been previously m easured using 
groun d-based high-precision spectroscopy (e.g., Beddi ng" et alj 
2010b) and the sp ace-based photometric mission CoRoT (e.g., 
Mic hel et alJl2008h . During the first seven months of Kepler sci- 
ence operations, an asteroseismic survey of solar-type stars made 
it possible to detect solar-like oscillations in about 500 targets 
(Chaplin etai1 l2011bl) . This constitutes an increase of one order 
of magnitude in the number of such stars with confirmed oscilla- 
tions. This large, homogeneous data sample opens the possibility 
of conducting ensemble asteroseismology. 

Since the start of Kepler science operations in 2009 May, 
a selection of solar-type stars have been continuously moni- 
tored at short cadence for more than seven months in order to 
test and validate the time-series photometry. Such continuous 
and long observations, also previously achieved by CoRoT (e.g., 
iBenomar et alj l2009bl) . are unprecedented for solar-type stars 
other than the Sun. We present herein the analysis of two of 
these starsfl namely, KIC 10273246 and KIC 10920273, both 
displaying relatively low SNR in the p-mode peaks. Two other 
stars, namely, KIC 11 395018 and KIC 1 1234888, are analysed 
in a companion paper dMathur et alJ201 ll) . The analysis of these 
four stars has been conducted in a way so as to group together 
stars observed for the same length of time. 

The two solar-type stars selected for this study are relatively 
faint (see Table [TJ if we bear in mind that the apparent mag- 
nitude target range for detection of solar-like os cillations with 
Keple r spans Kp as 6.5 to Kp sa 12.5 (see fig. 5 of Chaplin et al 



201 lal). The Kepler Inpu t Catalog (KIC; e .g., lLatham et al . 
20051: iBatalha et all 120101 grown et all 1201 lb . from which all 



KASC targets have been selected, classifies KIC 10273246 as 
an F-type star and KIC 10920273 as a G-type star (Table [B- The 
atmospheric parameters provided by the KIC - as derived from 
photometric observations acquired in the Sloan filters - do not 
have sufficient precision for asteroseismology. Although we can 
apply scaling relations to convert the KIC parameters of these 
targets into predicted seismic and non-resonant background pa- 
rameters, as well as fundamental stellar properties, caution is 
needed if use is to be made of these derived quantities. Tighter 



1 Within the Kepler Asteroseismic Science Consortium (KASC), 
KIC 10273246 is referred to as "Mulder" and KIC 10920273 as 
"Scully". 



constraints on T e s , log g and [Fe/H] will be obtained from spec- 
tra collected for these two targets with the FIES spectrograph at 
the Nordic Optical Telescope (Creevey et al., in preparation). 

The outline of the paper is as follows: We start in Sect. [2] 
by providing some background information on the properties of 
solar-like oscillations. In Sect. [3] we give an overview of the dif- 
ferent peak-baggingP| strategies employed and define a recipe for 
selecting frequency lists. Section|4]is devoted to a thorough anal- 
ysis of the power spectra of the time series. A summary and con- 
clusions are presented in Sect. [5] 

2. Properties of solar-like oscillations 

Solar-like oscillations are predominantly global standing acous- 
tic waves. These are p modes (pressure playing the role 
of the restoring force) and are characterized by being in- 
trinsically damped while simulta neously stochastically excite d 
by near-surface convection (e.g., Christensen-Dalsgaard 2004). 
Therefore, all stars cool enough to harbor an outer convective en- 
velope - whose locus in the H-R diagram approximately extends 
from the cool edge of the Cepheid instability strip and includes 
the red giant branch - may be expected to exhibit solar-like os- 
cillations. 

Modes of oscillation are characterized by three quantum 
numbers: n, I and m. The radial order n characterizes the be- 
havior of the mode in the radial direction. The degree I and the 
azimuthal order m determine the spherical harmonic describing 
the properties of the mode as a function of colatitude and lon- 
gitude. In the case of stellar observations, the associated whole- 
disk light integration and consequent lack of spatial resolution 
strongly suppress the signal from all but the modes of the low- 
est degree (with / < 3). For a spherically symmetric non-rotating 
star, mode frequencies depend only on n and I. 

The observed modes of oscillation are typically high-order 
acoustic modes. If interaction with a g-mode (gravity playing 
the role of the restoring force) can be neglected, linear, adia- 
batic, high-order acoustic modes, in a spherically symmetric star, 
satisfy an asymptot ic relation for the frequencies (Vandakurov 
l!967HTassoull [l980): 



v n i ~ Av (n + 1/2 + e)-/(Z+ l)D , 



(1) 



where the large frequency separation Av is the inverse of the 
sound travel time across the stellar diameter, s is a phase mostly 
sensitive to the properties of the near-surface region, and Do is 
a parameter sensitive to the sound-speed gradient near the core. 
The regular spacing of the frequency spectrum as conveyed by 
Eq. Q] is a characteristic feature of solar-like oscillations. We 
should, however, bear in mind that Eq.[T]is only an approxima- 
tion. The large frequency separation does in fact depend both on 
frequency and on mode degree, being defined as 



Av„( = v„, - v„_i;«Av. 



(2) 



It is also conventional to define a so-called small frequency sep- 
aration, also varying with frequency: 



Sv„i = v„i - y„_i i +2 ~ (4/ + 6)D . 



(3) 



The large frequency separation essenti ally scales with the squar e 
root of the mean stellar density (e.g.. lBrown & Gillilandl [l994). 
Furthermore, the small frequency separation is sensitive to the 
structure of the core, decreasing with increasing stellar age. 



3 The term peak-bagging refers to the extraction of individual mode 
http: //archive . stsci .edu/kepler/kepler_fov/search.php parameters from the power spectrum of a light curve. 
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Table 1. Information given in the Kepler Input Catalog. 



Star 


2MASS ID 


Kp 




Recalibrated" T e g 


logg 


[Fe/H] 


R 








(K) 


(K) 


(dex) 


(dex) 


(Rq) 


KIC 10273246 


19260576+4721300 


10.90 


6074 + 200 


6380 ± 76 


4.2 ± 0.5 


-0.3 ± 0.5 


1.506 


KIC 10920273 


19274576+4819454 


11.93 


5574 ± 200 


5880 ± 53 


4.1 ±0.5 


-0.4 ± 0.5 


1.594 



Notes. 

{a) A recalibration of the KIC photometry has posteriorly been performed by Pinsonneault et al. (in preparation). 



These two quantities thus have great diagnostic pote ntial (e.g., 
IChristensen-Dalsgaardlll993l:lDeheuvels et alj|2010bl) . 

Sharp variations in the stellar interior cause detectable os- 
cillatory signals in the frequencies, also visible in the behav- 
ior of Ay as a function of frequency ( e.g.. iMonteiro et alJ [2000: 
Ballot et all l2004t iBasu et all HQQ4; Cunha & Metcalfe 2007; 



Houdek & Gou gh 2007). These sharp features are mainly linked 



to borders of convection zones and to regions of rapid varia- 
tion in the sound speed due to ionization of a dominant ele- 
ment. Their combined signature is detectable in frequencies of 
low-degree modes and such an analysis becomes possible in the 
stellar case once frequency precision is sufficiently high. 

Stellar rotation, as well as any other physical process result- 
ing in departure from spherical symmetry, introduces a depen- 
dence of the frequencies of non-radial modes on m. When the 
cyclic rotational frequency of the star, v rot , is small and in the 
case of rigid-body rotation dominated by advection, t he cyclic 
frequ ency of a non-radial mode is given to first order by dLedouxl 
1951): 



Vnlm = VnlO + m V rot , M < I . 



(4) 



To a second order of approximation, centrifugal effects that dis- 
rupt the equilibrium structure of the star are taken into account 
through an additional frequency perturbation (independent of the 
sign of m). This perturbation in turn scales as the ratio of the 
centrifugal to the gravitational forces at the stellar surface, i.e., 
Q. 2 R 3 /(GM), where Q. denotes the surface angular velocity, R the 
radius of the star, M its mass, and G the universal gravitational 
constant. Although negligible in the Sun, these e ffects may b e 
significant for faster-rotating solar-type stars (e.g., Ballot 2010). 
Large-scale magnetic fields may also introduce further correc- 
tions to the oscillation frequencies. 

The frequency dependence of the mode surface amplitudes is 
determined both by (i) the frequency dependence of the stochas- 
tic process of excitation (mode energies result from a balance 
between the frequency-dependent energy input and the damp- 
ing rate) and by (ii) t he mode properties i n the region of vig - 
orous convection (e.g jHoudek et al.lll999l :fSamadi et al. 2007). 
The stochastic process of excitation is characterized by a rel- 
atively slow variation with frequency, meaning that it excites 
modes over a large frequency interval to comparable surface 
amplitudes. At low frequencies modes are evanescent in the 
region of efficient excitation, leading to small surface ampli- 
tudes. At high frequencies - greater than or equal to the acous- 
tic cut-off frequency - modes undergo considerable energy loss 
through running waves in the atmosphere. Excitation is most ef- 
ficient for those modes whose periods match the timescale of 
the near-surface convection. Also, the frequency of maximum 
amplitude, y max , is supposed to scale with the acoustic cut-of f 
frequency, v ac dBrown et al ] |!99lt iKieldse n & Bed ding|ll995l) . 
All this gives rise to a characteristic distribution of power with 
frequency which is a signature of the presence of solar-like os- 
cillations. 



Substantial changes in the properties of solar-like oscilla- 
tions occur with stellar evolution, particularly following the ex- 
haustion of hydrogen in the core. Most noticeable is the occur- 
rence of avoided cr ossings due t o coupling between p and g 
modes of like degree dOsakill975l :[Aizenma n et alj 19771) . which 
lead to significant departures from the regular frequency spacing 
described by Eq. Q] in the case of evolved stars. The frequen- 
cies of non-radial modes, in particular those of I = 1 modes, are 
shifted by avoided crossings when they couple with g modes 
trapped in the deep stellar interior. At the avoided crossings these 
modes have a mixed nature, with both p- and g-mode behav- 
ior. Provided they are excited to observable amplitudes (their 
high mode inertia reduces their surface amplitude), these so- 
called mixed modes are of great diagnostic potential because 
they probe the stellar core and are very sensitive to stellar age. 



3. On extracting estimates of mode frequencies 

3.1. Overview of the different fitting strategies 

We computed the power density spectrum (PDS) of the 
time series based on the implementation of the Lomb- 
Scargle periodogram (|Lomb1ll976t IScargldfl982l) presented in 
iPress & Rvbickil ([T989). This algorithm carries out reverse in- 
terpolation of the data onto a regular mesh and subsequently 
employs the fast Fourier transform. The power spectrum was 
then calibrated so that it satisfies Parseval's theorem, i.e., so that 
the total power in the positive-frequency side of the spectrum is 
equal to the variance of the time series (single-sided calibration). 
The effect of the window function is further taken into account 
when normalizing the PDS. 

A total of eleven individual fitters (A2Z.CR, A2Z_DS, 
A2Z_RG, AAU, IAS_OB, IAS_PG, IAS.TA, OCT, ORK, QML 
and SYD) extracted estimates of the p-mode frequencies for at 
least one of the two stars and subsequently uploaded their results 
to the Cat Basket data exchange facility. Different fitting strate- 
gies have been adopted and sometimes the same fitting strategy 
has been applied in an independent manner. All the fitting strate- 
gies adopted are, however, based on Fourier methods, the main 
idea behind them being either the maximization of the likeli- 
hood of a multi -parameter model describing the data or a classic 
prewhitening method. 

A frequency-domain representation of the data aims at mod- 
eling the limit PDS of the time series. Such a model typically in- 
cludes a sum of symmetric Lorentzian profiles meant to describe 
the individual p modes, together with a flat term and a number of 
additional ter ms describing both i nstrumental and stellar back- 
ground noise (lAnderson et alj|1990t) : 



6i 1 + P( v 



H, 



nlm 



[2(v - v„i m )/r nlm ] 



+ B(v) , 



(5) 



http : //bi son . ph . bham . ac . uk/kcatbaske t/ 
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where H is the mode height, F is the mode linewidth (related 
to the mode lifetime or amplitude e-folding time, r mo d e , through 
jtY = l/T mt) deX an d B(v) represents the background signal. The 
components arising from the decay of active regions, granula- 
tion an d faculae are commonly represent ed using a Harvey-like 
model dHarvevll985l:lAigrain et afll2004h : 



Zj 1+(2wT{)" 



+ W, 



(6) 



where {cr^} are the amplitudes, {t^} are the characteristic 
timescales, [sk] are the slopes of the individual power laws in 
the denominator, and W is a constant representing white noise 
(mainly due to photon shot noise). Equations [5] and [6] are rep- 
resentative of the models employed and individual fitters were 
further allowed to customize their own models. Hen ce we can- 
not tal k of a reference model, as was done in lAppourchaux et alj 
(120081) . 

Most fitters opted for performing a global fit (e.g., 
Appourc haux et alj|2008l) - whereby the whole set of free pa- 
rameters needed to describe the observed spectrum was opti- 
mized simultaneously - whil e there was still room for pseudo- 
global (or local) fitting (e.g., Jim enez-Reves et aT]|2008l) . an ap- 
proach traditionally adopted for Sun-as-a-star data, whereby nar- 
row frequency windows are considered at a time. 

Statistical inference from the data was either frequentist or 
Bayesian in nature. The former approach relies on the straight- 
forward application of a maximum likelihood estimator (MLE) 
taking into account the correct statistics of the spectrum (e.g., 
iTouta in & Appourc hauxll 1 9941 1 Appourchaux et al.ll2008l) or, in 
the special case of Gaussian statistics, a least squares estima- 
tor (LSE). To be specific: A2Z_CR applied a global LSE to 
a smoothed version of the spectrum; A2Z_DS applied a local 
MLE; IAS.TA, OCT and QML applied a global MLE. The lat- 
ter approach makes it possible to incorporate relevant prior in- 
formation through Bayes' theorem. A regularize d version of a 
MLE , known as maximum a posteriori (MAP; e.g.. Gaulm e et alJ 
2009), is a Bayesian point estimation method and was applied 
(globally) by A2Z_RG and IAS.PG. AAU and IAS_OB also per- 
formed global fits but instead used Markov chain Monte Carlo 
(MCMC) techniques to map the posteri or probability distribu- 
tions of the frequency parameters (e.g.. iBenomar et all [2009 a ; 
iGruberbauer et al.l2009tlHandberg & Campantel201 lh . With the 
exception of AAU and IAS_OB, who estimated the frequency 
uncertainties based on the posterior distributions of these same 
parameters, all the remaining fitters provided formal error bars 
derived from the inverse Hessian matrix. 

ORK and SYD adopted a different approach, which did 
not involve fitting a model to the power density spectrum. 
ORK applied a pr ocedure known as i t erative sine-wave fit- 
ting (ISWF; e.g. , iBedding et all l2007t iBonanno et aD 120081: 
White et a 02010h : it iteratively removes sinusoidal components 
from the data, which are identified as the maxima of the Fourier 
spectrum of the residuals. Since the number of oscillation modes 
present in the data is unknown a priori, the algorithm requires a 
stopping rule that associates some degree of confidence to the 
amplitude of each extracted sinusoidal component. For the anal- 
ysis by SYD, the power spectrum was smoothed by convolving 
with a Gaussian with full width at half maximum of 1 .4 fiHz 
(comparable to the intrinsic linewidth of the modes) and the 
frequencies of the highest peaks were measured. Mode identi- 
fication (assigning I values) was done using the echelle diagram 
dGrec et al.iri983h 'both by ORK and SYD. Finally, frequency 



uncertainties are estimated either analytically (ORK) or by per- 
forming Monte Carlo simulations with artificial data (SYD). 

3.2. Procedure for selecting frequency lists 

A procedure for selecting relevant frequ ency lists to be used i n 
stellar modeling has been described in Metcalfe et al. (2010). 
However, we have opted for revising that procedure. This has 
been motivated, at first, by the fact that different fitting strate- 
gies provided frequency uncertainties that differed greatly. For 
instance, uncertainties that had been underestimated would un- 
desirably upweight the respective individual frequency sets. 
Consequently, a more robust method of frequency selection was 
needed. Other refinements were also made and we give below a 
detailed description of the new, revised procedure. 

We aim at selecting, for a given star and for a given mode 
degree identification, two frequency sets - a minimal frequency 
set and a maximal frequency set - that will provide initial con- 
straints for the modeling and allow for refined model-fitting, re- 
spectively. As before, we would like to determine individual sets, 
as opposed to averaged sets, meaning that these sets are fully re- 
producible. 

We start by constructing both a minimal list and a maximal 
list of modes. Take N as being the total number of individual 
fitters providing peak-bagging results for a given star and as- 
sume for now that N >2. For each {n, /} pair, we a pply Peirce's 
criterion (see Appendix lAl for its implementation: |Peircdll852t 
Gould 1855) for the rejection of outliers and assess how many 
frequency estimates are retained. Inclusion or not of the mode 
in the minimal list then results from a vote including all (N) fit- 
ters: If the number of frequencies retained is greater thai£| [N/2} 
then the mode is added to the minimal list. Inclusion or not of 
the mode in the maximal list results instead from a vote only 
including fitters providing estimates for that particular mode: If 
the number of frequencies retained is at least 2 then the mode 
is added to the maximal list. The minimal list is thus a subset of 
the maximal list. For N — 2 the minimal and maximal lists are 
degenerate and will coincide. 

In the final stage of the procedure, we compute for each 
of the N individual frequency sets the normalized root-mean- 
square deviation (nrmsd) with respect to the frequencies aver- 
aged over all contributing fitters, {v„/}, belonging to the minimal 
list of modes: 



nrmsd 



ZO„/ - v«;) 2 /o"„, 

tl,l 



(7) 



where cr„i is the uncertainty in v„i, and Ni is the number of 
modes in a particular individual set that actually belong to the 
minimal list. The best fit is defined as being the individual set 
with the smallest normalized rms deviation. Note that by best 
fit we mean here the most representative fit among the N avail- 
able sets and not necessarily the one closest to the truth. The 
minimal and maximal frequency sets are finally given by those 
modes provided by the best fit that belong to the minimal and 
maximal lists, respectively. The homogeneous character of the 
revised procedure is reassured by generating the minimal and 
maximal frequency sets from the same individual set. The cor- 
responding frequency uncertainties are simply those associated 
with the best fit. When N = 2 a single frequency set is defined, 
the same happening for N- 1 in which case it coincides with the 
only available individual set. 



LxJ returns the closest integer less than or equal to x. 
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Although having been revised, this procedure is still subject 
to future improvement. The main reason for this is the fact that 
we assume the I values to be correctly assigned to modes by all 
the fitters. A way of overcoming this would be to implement a 
clustering algorithm that groups frequencies and posteriorly tags 
them with a mode degree based on a vote. Another issue con- 
cerns the possibility that the best fit misses some modes from 
the minimal and maximal lists. A solution for this would be to 
ask the fitter providing the best fit to reanalyse the power spec- 
trum. One last drawback of this procedure is the possibility that 
the minimal and maximal frequency sets may contain some out- 
liers, i.e., frequency estimates that have been ruled out according 
to Peirce's criterion. 



4. Data analysis 

4.1. Time series preparation 

The stars were observed from 2009 May 2 to 2010 March 19, 
i.e., from Quarter (Q0) to Quarter 4 (Q4). The duty cycle 
over the course of the approximately 10 months of observations 
was above 90%. Available time series suffered from several in- 
strumental perturbations , and we have thus decided to develop 
our own corrections (Garcf a~et al.l 1201 it) , independently from 
the pre-s earch data condition ing module of the Kepler Science 
Pipeline (Jenk ins et al.ll2010l) . which generates corrected light 
curves for transiting-planet search. 

The raw-flux time series were corrected for three types of ef- 
fects: outliers, jumps and drifts. Data points considered to be out- 
liers exhibit point-to-point deviations greater than 30"^, where 
(Tdift is the standard deviation of the first differences of the time 
series. This correction removed approximately 1% of the data 
points. Jumps are sudden changes in the mean value of the time 
series caused by, e.g., attitude adjustments or a drop in pixel sen- 
sitivity (see Fig. [TJ. They were identified as causing spurious 
differences in the mean power of contiguous bins spanning one 
day. In cases where we know the photometric apertures to have 
changed, we by default check for jumps (only implemented in 
segments Q2.1 and Q4.1). Each jump has been manually vali- 
dated. Finally, drifts are small, low-frequency perturbations due 
to temperature changes (e.g., after a long safe-mode event) and 
lasting for a few days (see Fig.|2]i. The corrections are based on 
the software developed to deal with high- voltage p erturbations in 
the GOLF/SoHO instrument dGarcfa et al.ll2005l) . We fit a 2nd- 
or 3rd-order polynomial function to the section of the time series 
where a thermal drift has been observed after comparing several 
light curves of the same Quarter. The fitted polynomial is then 
subtracted and another polynomial function of 1st- or 2nd-order 
- used as a reference - is added, which has been computed based 
on the observations done before and after the affected section. If 
the correction has to be applied to a border of the time series, 
then only one side of the light curve is processed. 

Once the aforementioned corrections have been applied, we 
merge the data of the different Quarters into a single time se- 
ries, after equalizing the average counting-rate level between the 
Quarters (or sometimes even within a Quarter when some instru- 
mental parameters have been changed). In order to do so, and to 
allow conversion into units of ppm, we use a series of 6th-order 
polynomial fits, one for each segment. Finally, we normalize the 
standard deviation of the data obtained during Q3 to the aver- 
age of the other Quarters, since Q3 is considerably noisier for 
the two stars. This normalization proved to be a good compro- 
mise between using the noisy Q3 data and not using them at all. 
Given the insufficient technical information available, it is not 




Fig. 1. Raw-flux time series corrected only for outliers (top 
panel) and corrected relative time series (bottom panel) for the 
case of KIC 10273246. Vertical dashed lines mark the beginning 
of each Quarter, whereas vertical dotted lines separate the seg- 
ments within a Quarter. 




80 

Time (days 



Fig. 2. Raw-flux time series (corrected only for outliers) of 
Quarter 2 for the case of KIC 10273246. Vertical dotted lines 
separate the three segments of Q2. The induced thermal drift af- 
ter a safe-mode event is visible between the 60th and 70th days. 



possible to know whether or not there is a variation in the gain 
of the CCD module during Q3. If such a gain variation is indeed 
present, the mode amplitudes would be affected and the adopted 
normalization will then correct these amplitudes toward some- 
thing closer to their real value. If, on the other hand, no gain 
variation is present, then the adopted correction will effectively 
reduce the mode amplitudes. Having reduced the standard de- 
viation of the data in Q3 by 34% for KIC 10273246 and 49% 
for KIC 10920273, and given that Q3 represents 3/7 of the total 
length of the time series used for peak-bagging (i.e., from Q0 to 
Q3; see Sect. 14.51 ). then mode amplitudes would be reduced by 
about 14% and 21%, respectively. However, we believe this lat- 
ter scenario to be less likely and hence the quoted values of 14% 
and 2 1 % can only be regarded as upper limits to mode amplitude 
reduction, which in reality may be much smaller. 
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Fig. 3. Low-frequency end - between 0.1 and 3 //Hz - of the Fig. 4. Similar to Fig. [3]but for the case of KIC 10920273. 
power density spectrum of KIC 10273246. 



4.2. Rotational modulation 

To investigate the stellar rotational period (P m t) of these stars 
we seek high-SNR peaks in the low-frequency end of the PDS. 
However, the procedure for merging the different data sets de- 
scribed in the previous paragraph filters out the power density 
below 1 //Hz. Therefore, starting with the corrected data, we have 
generated new merged time series in which a triangular smooth- 
ing - over a period selected from the range of 12-20 days - has 
been used to normalize the ligh t curves instead o f the 6th-order 
polynomial fit (for details, see iGarcfa eta! ] 120111) . Finally, we 
computed the PDS of the full-length time series and also of two 
subsets of 160 days each. The low-frequency ends of the PDS are 
shown in Figs. |3j and g] for KIC 10273246 and KIC 10920273, 
respectively. A different pattern of peaks appears for each star. 
In the case of KIC 10273246 (Fig. [3), the highest peak is seen at 
about 0.50 //Hz (P mt as 23 days) in the three spectra considered, 
without any signature of differential rotation. This turns out to 
be a reliable signature of rotational modulation, since compari- 
son with the power spectra of other stars observed with the same 
CCD module ruled out the possibility of an instrumental arti- 
fact. On the other hand, KIC 10920273 (Fig. g} exhibits a high 
peak at about 0.43 //Hz, though mostly during the second subset 
of 160 days. Indeed, this peak is not very significant, and during 
the first subset it is at the same level as that of the adjacent peaks. 
The PDS of the full-length time series also reveals itself as being 
rather noisy, and it is difficult to disentangle the stellar rotational 
period from it. Therefore, while a cyclic rotational frequency of 
0.43 //Hz (P Iot « 27 days) seems plausible, we cannot exclude 
the possibility that this peak is associated with the roll sched- 
ule of the telescope. One of several alternative explanations for 
the low-SNR peaks would be that this star was observed during a 
low magnetic-activity period, with a small number of spots on its 
surface. A deeper study based on a longer data set is necessary 
to confirm the rotational rate of KIC 10920273. 



4.3. Power spectral density of the background 

The power density spectra - based on the first seven months 
of observations, i.e., from Q0 to Q3 - of KIC 10273246 and 
KIC 10920273 are displayed in Fig. [5] The background signal 
rises toward the low-frequency end of the spectra with contribu- 
tions from granulation and activity. At the high-frequency end, 



the spectra are dominated by photon shot noise. Also visible at 
the high-frequency end are a number of peaks corresponding 
to harmonics of the inverse of the long-cadence period, an ar- 
tifact appearing in pow er spectra of short-cadence time series 
( Gilli land et al.l 1201 Obi) but having negligible influence on the 
current study. In between these two frequency regimes there is 
a conspicuous cluster of power due to the presence of p modes, 
with the hotter and brighter target KIC 10273246 displaying the 
higher SNR. 

We fitted a model similar to the one described in Eq. [6] 
to a heavily smoothed version of both spectra. This model in- 
cluded an additional Gaussian fun ction aimed at describin g the 
p-mode power-excess hump (e.g., Kalling er et alJl20 10b). The 
fitting window was from 100 //Hz up to the Nyquist frequency, 
hence we did not consider a term accounting for the decay of 
active regions, whose typical timescale is considerably longer. 
We found no evidence for a facular component based on a sim- 
ple visual inspection of the spectra, thus having not included 
such a component in the background model. Faculae had pre- 
vious ly been reported o n a couple of Kepler solar-type targets 
bylChaplin et al] d2010l) . A component carrying the signature of 
stellar granulation is, however, clearly displayed by both stars. 
Figure[5]displays the fitted models, while values of the fitted pa- 
rameters are given in Table [2] Note that the slope of both granu- 
lation components has been fi xed at the same valu e, i.e., s gran = 3, 
which is closer to solar (e.g., Mi chel et al 1 12009) than the value 
of 2 originally proposed bv lHarvevi (11985 ). Although we lack a 



physical reason for fixing this slope, the intention is to facilitate 
the comparison between the values of Tg ran . 

The high-frequency noise power spectral density (W) is ap- 
proximately 14% and 16% higher than predictions in the cases of 
KIC 10273246 and KIC 10920273, respectively. These predic- 
tions have been computed according t o the empirical minima l 
term model for the noise presented in Gillilan d et al.l (|2010b), 
which takes into account the Kepler-band magnitude of the star 
and the performance of the instrument. Noise levels are thus 
close to being Poisson-limited. A likely source of this extra noise 
is the larger scatter in Q3. 



4.4. The global asteroseismic parameters Av, 



and Sv, 



nO 



In the past few years, a number of automated pipelines have been 
developed to measure global asteroseismic parameters of solar- 
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Fig. 5. Power density spectra - smoothed over 1 //Hz in order to enhance p-mode visibility - of KIC 10273246 (top panel) and 
KIC 10920273 (bottom panel), plotted on a log-log scale. The solid black lines represent the fits described in the text. The remaining 
lines represent different components of the fitting model: the power envelope due to oscillations added to the offset from white noise 
(dashed), and granulation (dotted). 



Table 2. Background model parameters describing non-resonant 
features. 



Star 


fgran 


fgran 


., a 


W 




(ppm) 


(s) 




(pprrr/iHz -1 ) 


KIC 10273246 


69.9 ± 0.8 


390 ± 12 


3 


4.15 ±0.02 


KIC 10920273 


80.1 ±0.4 


351 ±4 


3 


11.58 ±0.02 



Notes. 

^ The slope of both Harvey-like terms has been fixed. 



like o scillators dHuber et"al 
12009 

2OT0t lKaroffetaLM20T 



2009; iRoxburghl 120091: ICampante et al.1 l2010t iHekkeretai 



Table [3] we give a representative set of estimates of Av and 
Vmax for the two stars in our study. These estimates were sup- 
plied by one of the pipelines , namely, the one described in 
Mosser & App ourchauxl (120091) . and are based on the analysis 
of the full-length time series. The large separation is a function 
of frequency. Consequently, the extent of variation of its mean 
value will depend upon the variation in frequency range adopted 
for its computation. However, as long as this frequency range 
includes v max , the impact of small differences in the range can 
be neglected. For the sake of completeness, we give the adopted 
frequency ranges in Table [3] An observed relation between Av 
20091: Mosser & Appourchauxl and v max for solar- like oscilla tions in main-sequence stars is pre- 
sented in lStello et all d2009ah : 



Mathur et alj 



2010b) 



Most of them 

have already been successfully tested on CoRoT data (e.g., 
Math ur et al.ll2010"ah . The automated nature of these pipelines 
is required if we are to efficiently exploit the plenitude of data 
made available by Kepler on these targets. A thorough compar- 
ison of complementary analysis methods used to extract global 
asteroseismic parameters of main-sequence and subgiant solar- 
like oscillators is presented in Vern er et"a"fl (1201 ll) . 

The results of the different pipelines on global parameter 
extraction have also been uploaded to the Cat Basket. For a 
given global parameter, some groups submitted results using 
more than one method or on more than one data release. In 



Av oc V 



,0.77 



(8) 



The values quoted in Table [3] for these two global parameters 
satisfy this relation, as can be seen in Fig. [6] Also, T yran is seen 
to scale inversely with v max , as predicted by Kieldsen & Bedding 
(I20T1 . 

The same pipeline has provided estimates of 5v„o based on 
the analysis of seven months of data. These are also given in 
Table[3] 
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Table 3. Estimates of the global asteroseismic parameters Av, 
v max and 6v„ . 



Star 


<Av> 


Range 


"max 


(8v„ ) 




C"Hz) 


0"Hz) 


(pHz) 


OuHz) 


KIC 10273246 


48.2 ± 0.5 


[537, 1140] 


839 ±51 


5.6 ± 1.2 


KIC 10920273 


57.3 ±0.8 


[757, 1290] 


1024 ± 64 


6.0 ± 1.5 




2000 3000 
"„» C/^Hz) 



Fig. 6. Observ ed relation between Av and v max . Here we repro- 
duce fig. 2 of IStello et all (l2009al) in linear scale, where only 
stars whose v max is greater than the Nyquist frequency for long- 
cadence sampling are displayed. The solid line is the power-law 
fit in their eq. 1 with the corresponding ±15% deviations shown 
as dotted lines. We added to the plot the two stars being consid- 
ered in this s t udy (f ive-pointed stars), the two stars (squares) of 
Math ur et al. (120 ll). an d the three bright G-type stars (triangles) 
of lChaplin et all d2010t) . 



4.5. Mode frequencies 

Estimates of mode frequencies were obtained by several individ- 
ual fitters for both stars, based on the first seven months of obser- 
vations (from Q0 to Q3). The procedure described in Sect. 13.21 
was then used to select minimal and maximal frequency sets. 
An ad hoc step was, however, implemented at this stage, that 
aimed at removing doubtful modes from these sets based on 
SNR considerations. This quality control led to the removal of a 
few modes at the low- and high-frequency ends of the sets. 

Another issue relates to the presence of mixed modes in the 
power spectrum. As we will see below, the best fit relies, for both 
stars, on an algorithm based on a frequency-domain representa- 
tion of the data (cf. Eq. [5]). Moreover, the deterministic models 
that have been used only contained modes of degree up to / = 2, 
with three modes per radial order, one of each degree. The pres- 
ence of an extra mixed mode in the vicinity of a mode with low 
SNR therefore greatly complicates the fit due to the inadequacy 
of the model. As a result, we opted for prewhitening the longer- 
lived mixed modes and give their frequencies a posteriori, while 
discarding any neighboring power structure with low SNR fitted 
by the MLE-based algorithm. 

Eleven individual fitters (A2Z_CR, A2Z_DS, A2ZJRG, 
AAU, IAS_OB, IAS.PG, IAS_TA, OCT, ORK, QML and SYD) 
have provided results for KIC 10273246. A summary of the 
anonymous best fit is given in Table [4] The fitter has also im- 
posed priors on the parameters used to describe the background 
signal, which were based on the output of a previous fit to the 



background alone. Tableland Fig. [7] combine the minimal and 
maximal frequency sets (recall that the former is a subset of the 
latter). We have identified a total of 30 modes of degree / = 0, 1 , 2 
that span a frequency range of width * 10 Av. The mode at 
1122.70 //Hz has only been reported by two fitters, who have 
identified it as being a radial mode despite its alignment with the 
1 = 2 ridge. Given its very low power and the known convergence 
issues experienced by MLE-based methods under such low SNR 
conditions, we have decided to change its tagging a posteriori 
from I = to I = 2. This mode should thus be considered with 
caution. 

The problem of ridge identification (i.e., the tagging of 
modes by degree /) in F stars dates back to when CoRoT obser- 
vation s of HD 49933 were first analysed by Appourchaux et all 
(2008j). In the present case such an identification can be sim- 
ply done by visual inspection of the echelle diagram in Fig. [7] 
Nonetheless, two of the fitters (IAS.OB and IAS.TA) provided 
results for the two complementary identification scenarios. From 
the results provid ed by IAS-TA we are able to compute the like- 
lihood ratio (e.g..lAppourcha ux et al.l lT998) in favor of our pre- 
ferred scenario, whereas computation of its Bay esian counter- 
part, the more conservative Bayes' factor (e.g.. iLiddlel [2009: 
lHandberg & Campantdl20i"il) . is possible based on the results 
provided by IAS_OB. Both approaches returned conclusive val- 
ues in support of the adopted scenario. 

Ten individual fitters (A2Z_CR, A2Z_DS, A2Z_RG, 
IAS.OB, IAS_PG, IAS.TA, OCT, ORK, QML and SYD) 
have provided results for KIC 10920273. A summary of the 
anonymous best fit is given in Table [4] In addition to the 
information provided in that table, it should also be mentioned 
that a Gaussian prior has been imposed on 6v„q. Tableland 
Fig. [8] combine the minimal and maximal frequency sets. We 
have identified a total of 21 modes - considerably less than 
for KIC 10273246 - of degree I = 0, 1, 2 that span a frequency 
range of width m 8 Av. The / = 1 mode at 1135.36 //Hz is 
very close to the second harmonic (2 x 566.391 /iHz) of the 
long-cadence sampling. Although this anomaly is known to be 
more prominent for intermedi ate harmonics, we no netheless 
follow the recommendation of Gillil and et al.l (1201 Obi) to flag 
this mode as suspect. 

The quasi-regularity of the small (<5v„o) and large frequency 
separations (Av„o and Av„2) is evident from Figs.|7]and|8] Notice 
that if these stars were to strictly obey the asymptotic rela- 
tion in Eq. [T] they would then exhibit vertical ridges in the 
echelle diagram provided use of the correct Av. The small sepa- 
ration <5v„o is, however, more clearly distinguished in the case of 
KIC 10920273, which might be an indication of smaller mode 
linewidths in this cooler star (see Sect. 14.61 for a discussion on 
mode linewidths). 

A striking feature in both echelle diagrams is the jagged 
appearance of the I = 1 ridge, a trademark of the presence 
of avoided crossings and an indicator of the evolved nature of 
these stars. These same features have also be en seen in the 
cases of ground -based observations o f tj Boo dKieldsen et al.l 
120031). B Hvi dBedding et al.l 120071) and possibly Procyon 
(Bedding et al. 2010b), as well as in the cases of the CoRoT tar- 
get HD 49385 dDeheuvels et al.ll2010al), a nd KASC survey tar- 
gets KIC 1 1026 764 dMetcalfe et al.l|2010h . KIC 11395018 and 
KIC 11234888 dMathur et all 1201 1 ). Figure M disp lays a so- 
called p-g diagram as introduced by Bedding (201 1), where the 
frequencies of the avoided crossings (i.e., the frequencies of the 
pure g modes in the core cavity) for a number of stars are plotted 
against the large separation of the p modes. Much of the diag- 
nostic potential of mixed modes can be captured in this way, 
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Table 4. Summary of the best fits. 



Star 


Method 


Height (H) 


Linewidth (r) 


Splitting (v rol ) 


Inclination" (0 


KIC 10273246 


Global 
MAP 


One free parameter (H/ = o) per order 
No prior 
Hi =l /H l=0 = 1.5 (fixed) 
Hi =2 /H, =0 = 0.5 (fixed) 


One free parameter per order 
No prior 


Free 
No prior 


Free 
No prior 


KIC 10920273 


Global 
MAP 


One free parameter (H t=0 ) per order 
No prior 
H l=l /H l=0 = 1.5 (fixed) 
Hi =2 /Hi =0 = 0.5 (fixed) 


One free parameter per order 
Gaussian prior 


Fixed at /jHz 


Fixed at 0° 



Notes. 

(a) Introduced in Eq.fBI 
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t 300 
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Frequency modulo 48.50 {jlHz) 

Fig. 7. Echelle diagram of the power density spectrum of 
KIC 10273246 (the colorbar decodes the power density level). 
The minimal (filled symbols) and maximal (filled and open sym- 
bols) frequency sets are displayed. Symbol shapes indicate mode 
degree: 1 = (circles), /= 1 (triangles) and 1 = 2 (squares). 



since their overall pattern is determined by the mode bumping 
at each avoided crossing, which in turn is determined by the g 
modes trapped in the core. This diagram could prove to be an 
instructive way to display results of many stars and to allow for 
a first comparison with theoretical models. We also report here 
the possible presence of a I = 2 mixed mode in the power spec- 
trum of KIC 10920273 (at 873.10jifHz) that should, however, be 
confirmed by stellar models. 

Detection of / = 3 modes with photometric observations 
is made very difficult due to geometric cancellation effects. 
Solar-like oscillations with / = 3 from Kepler photometry have 
nonethe less been repo r ted for a set of low-luminosity r ed gi- 
ants by iBedding et all d2010al) . iDeheuvels etal] d2010al) have 
also reported the presence of I — 3 modes for the CoRoT tar- 
get HD 49385. We should bear in mind that, except for ORK 
and SYD, all the remaining fitters used deterministic models in 
their frequency-domain representations of the data that only con- 
tained modes of degree up to I — 2, meaning that a statistical 
assessment of the presence or not of / = 3 modes could not be 
done. ORK and SYD, which were the only fitters that did not 
make any prior assumptions about the degree of the modes, have 
not reported the detection of modes that could be interpreted as 
1 = 3 modes. 



Table 5. The minimal and maximal frequency sets of observed 
oscillation frequencies for KIC 10273246. 



/ 


Frequency 


Uncertainty 




(/"Hz) 


O^Hz) 





737.90 


0.30 





785.40 


0.20 





833.90 


0.20 





883.50 


0.20 





932.70 


0.50 





981.10 


0.30 





1030.70 


0.40 





1079.30 


0.20 




622.80 


0.20 




661.90 


0.50 




695.75* 


0.27 




724.70 


0.20 




764.30 


0.30 




809.80 


0.20 




857.30 


0.20 




905.60 


0.30 




950.00 


0.30 




1008.60 


0.40 




1056.30 


0.20 




1103.30 


0.40 


2 


688.50 


0.70 


2 


734.80 


0.60 


2 


779.50 


0.40 


2 


830.30 


0.40 


2 


880.60 


0.50 


2 


927.50 


0.40 


2 


977.60 


0.40 


2 


1025.30 


1.30 


2 


1073.70 


0.20 


2 


1122.70 flr 


0.40 



Notes. 

<n) Mode belonging exclusively to the maximal frequency set. 

(b> 1=1 mixed mode introduced a posteriori. 

w Tagging changed a posteriori from 1 = to 1 = 2. 



4.6. A word on mode linewidths, heights, and amplitudes 

A thorough discussion of the mode linewidths, heights, and am- 
plitudes goes beyond the scope of this work. However, there are 
some aspects we would like to mention here. 

The intrinsic frequency resolution of the spectra (a;0.05//Hz) 
makes it possible to resolve the modes. This condi tion is obeyed 
provided the observation length T » 2r mot j e (Chaplin et al. 
2003). Figure [10] displays, for each star, the linewidths of the 
radial modes returned by the respective best fit (see also Tables 
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Table 6. The minimal and maximal frequency sets of observed 
oscillation frequencies for KIC 10920273. 



/ Frequency 
C"Hz) 



Uncertainty 
0"Hz) 



826.66° 
882.77 
939.58 
997.14 
1054.33 
1111.51 
1170.77" 
1226.34" 



0.25 
0.20 
0.16 
0.18 
0.30 
0.25 
0.33 
0.33 



794.65" 
838.61* 
914.52 
968.19 
1023.58 
1079.10 
1135.36"* 



0.32 
0.25 
0.16 
0.13 
0.14 
0.31 
0.31 



822.39" 
873.10" c 
934.49 
992.44 
1049.36 
1106.76 



0.28 
0.32 
0.22 
0.13 
0.39 
0.34 



Notes. 

(a) Mode belonging exclusively to the maximal frequency set. 
{b) I = 1 mixed mode introduced a posteriori. 
<c> Possible 1 = 2 mixed mode introduced a posteriori. 
(</) Mode very close to the second harmonic of the inverse of the Ion 
cadence period. 




20 30 40 50 

Frequency modulo 57.30 (fiHz) 



Fig. 8. Similar to Fig.|7]but for the case of KIC 10920273. 



[7] and [8). The radial modes considered are those belonging to 
the maximal frequency set. Notice the near-constancy with fre- 
quency of the mode linewidths in the case of KIC 10920273, 
whereas for KIC 10273246 the linewidths increase steadily, al- 
though with considerably larger error bars (note that both sets 
of error bars were derived from the inverse Hessian matrix). 
Using pulsation computations of a grid of stellar models and the 
first asterosei s mic re sults on mode lifetimes of solar-like stars, 
Chapl in et al .] (120091) suggested a simple scaling relation be- 



KIC 10273246 



KIC 11395018 



KIC 1 0920273 
KIC 11026764 HD 49385 



50 
Al/ C"Hz) 



Fig. 9. p-g diagram showing the frequencies of avoided crossings 
(diamonds) plotted against the large separation for a number of 
stars. The vertical lines show the range over which the / = 1 ridge 
is clearly seen. The inclined line is the observed relation between 
Av and v max of lStello et aD (l2009al) . For KIC 10273246 there are 
two clear avoided crossings, whereas for KIC 10920273 there is 
only one clear crossing plus another likely one (open symbol). 



tween the mean mode linewidth of the most prominent p modes 
and r eff : 



«- <r><xr e 4 ff . 



(9) 



Figure [TO] also displays the resulting predictions of mean 
linewidths of the most prominent modes (we have considered 
(T)q * 1.2 fjHz and r eff Q = 5777 K, and taken the recalibrated 
KIC temperatures). The agreement with the observed values is 
fairly good in the case of KIC 10273246. On the other hand, the 
predicted value obtained for KIC 10920273 using Eq.|9]overes- 
timates the observed linewidths. Assuming validity of this equa- 
tion, this might be the result of the combination of two factors, 
namely, an overestimation of and an overestimation of the 
fitted background that leads to underestimated fitted linewidths. 
Nonetheless, the F-type star KIC 10273246 exhibits the larger 
mode linewidths, as exp ected from Eq. [9] We should also note 
that lBaudin et al 1 ) found a much stronger dependence of 
the mean mode linewidth on T e ff . 

When a mode is resolved, as is the case here, it is the 
mode height, H, that determines the SNR in power, viz., the 
height-to-background ratio. Figure[TT]displays, for each star, the 
heights of the radial modes returned by the respective best fit 
(see also Tables [7] and [8). The modes are the same as shown in 
Fig. [10] We also indicate in Fig. [TJJ the SNR of the strongest 
radial modes, as well as the background fits of Fig. [5] (which 
are an indicator of SNR = 1). Overall, KIC 10273246 exhibits a 
higher SNR if we take into account the whole plotted frequency 
bands. Also apparent is the larger width of the p-mo de hump of 
KIC 10920273, which rou ghly scales with v max (e.g JStello et al.l 
120071; iMosser et al.ll2010t) . 

Furthermore, we compute d the rms amplitudes of the radial 
modes, A„o, according to (e.g., Chapli n et al . 2003}: 



- \l ■ 



(10) 



Amplitudes are always better constrained than the heights and 
the linewidths themselves. We were careful enough to com- 
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iKield sen & Beddind (1 19951) have suggested an empirical 
scaling relation to predict the amplitudes of solar-like oscil- 
lations that, although extensively used, predicts amplitudes 
in F-type stars th at are higher than actually observed (e.g., 
Mich el et alj |2008). Recently, the same authors have proposed 
a new scaling relation fo r the amplitudes whic h is based on sim- 
ple physical arguments (Kjeldsen & Bedding 20jjl): 



Lt 0.5 
mode 

eti 



(ii) 



where L is the stellar luminosity, and the value of r is chosen to 
be either r — 1 .5 (assuming adiabat ic oscillations) or r = 2 (fol - 
lowing a fit to observational data in Kj eldsen & Beddinglll995l) . 
By assuming that v max is a fixed fraction of the acoustic cut-off 
frequency, i.e., 



Fig. 10. Linewidths of the radial modes returned for 
KIC 10273246 and KIC 10920273 by the respective best 
fit. Modes represented by open symbols belong exclusively to 
the maximal frequency set. Dot-dashed lines mark the predicted 
mean linewidths of the most prominent modes using Eq. [9] 
The horizontal dimension of the line-filled areas represents the 
uncertainty in v max . 



oc v a , 



MT\ 5 

eft 



(12) 



and adopting a scalin g relation for the stellar mass based on seis- 
mic parameters (e.g., Kallinger et al. 2010b), 



(13) 




KIC 10273246 
KIC 10920273 



900 1000 
Frequency (/iHz) 



Fig. 11. Heights of the radial modes returned for KIC 10273246 
and KIC 10920273 by the respective best fit, plotted on a log- 
linear scale. Modes represented by open symbols belong exclu- 
sively to the maximal frequency set. Also indicated is the SNR 
of the strongest modes as well as the background fits of Fig. [5] 
(dotted lines). 



pute errors on the amplitudes that took into account the cor- 
relations between the fitted parameters F„o and H„o (e.g., 
App ourchauxl201 ll) . These amplitudes were then scaled to their 
bolometric equivalent using the bolometric correcti on derived 
from the spectral response of the Kepler passband dBallot et alJ 
1201 lb . Finally, we obtained a maximum bolometric amplitude 

of A 1om = 1X1 ± °- 75 PP m ( at 833 - 90 y« Hz ) for KJC 10273246 
(cf. TableE), and of A^j =7.56 ±0.99 ppm (at 997. 14 //Hz) for 
KIC 10920273 (cf. Table|8]l. These values were computed based 
on the results returned by the respective best fit and were found 
to be consistent with the values obtained from the other fitters' 
results. 



we can combine Eqs. [TjU - U3l to obtain the following relation (nor- 
malized with respect to values in the Sun): 



.(max) _ .(max) 
^nO.bol ,,0,bolO 



' eft' 



' eft'O 



1.5- 



-2.5 



Av 



Av ( 



_T_ 

r©, 



-0.5 



,(14) 



where ^0 



= 2.53 ppm, v max = 3050 yt/Hz, and Av0 = 
135 fiHz. Note that by setting r — 1.5 (as we will be assum- 
ing hereafter), the dependence of Eq.[14]on T e g is canceled and 
this scaling relation then solely depends on seismic parameters. 
Finally, use of this relation gives a predicted maximum bolomet- 
ric amplitude of A^, =7.27 ± 1 .33 ppm for KIC 10273246, and 

of A Sm = 8 - 24 ± lM PP m for KIC 10920273, with the rela- 
tively large uncertainties dominated by the errors on v max . These 
values agree with the observed values at the 1-cr level. This is 
a particularly interesting result, since no discrepancy is seen be- 
tween the predicted and observed amplitudes for the F-type star 
KIC 10273246. 



4.7. Variation of Av with frequency 

An independent methodology has been used to estimate the large 
frequency separation and the mode frequencies, which is based 
on the analysis of the spectrum with the envelope a utocorre- 
lation function (EACF; Mosser & Arjrjourchaux 2009). As ini- 
tially proposed by iRoxburgh & Vorontsovl d2006l) . the autocor- 
relation of the time series - or, equivalently, the power spectrum 
of the power spectrum - windowed with a narrow filter giv es the 
variat ion of the large separation with frequency, Av(v). [Mosser 
(2010) has shown that, with a dedicated comb filter for analyz- 
ing the power spectrum, it is possible to obtain independently 
the values of the large separation for the odd (/ = 1) and even 
(I = 0, 2) ridges. We denote them by Av ot jd and Av even . Proxies 
of the mode frequencies can then be integrated from the Av(v) 
frequency pattern. In practice, they are derived from the correla- 
tion between the observed PDS and a synthetic spectrum based 
on the Av(v) pattern. 
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Fig. 12. Variation of Av dd and Av even with frequency for 
KIC 10273246. The dashed line corresponds to the odd ridge, 
whereas the solid line corresponds to the even ridge. Error bars 
are given with the same linestyles. Proxies of the mode frequen- 
cies derived with the EACF have been superimposed and their 
values appended: / = (circles) and I = 1 (triangles). The inset 
shows the comb filter used in the analysis. 

The values of Av dd and Av even for KIC 10273246 are given 
in Fig.Q~2] Notice the agreement between the proxies of the fre- 
quencies derived with the EACF and the frequencies of the max- 
imal frequency set (Table[5]l. The EACF emphasizes the low val- 
ues of Av dd as well as its large gradient at the low-frequency 
end. Also clear is the large discrepancy relative to a regular 
echelle spectrum around 1000 juHz. The only mode present in 
the maximal frequency set and not detected with the EACF is 
the mixed mode at 695.75 frtiz. This peak appears as supernu- 
merary when compared to the regular agency of the modes. The 
EACF makes it possible to derive the large separation one radial 
order further than does peak-bagging. 

Results for KIC 10920273 are given in Fig. Qj] The lower 
SNR is counterbalanced by using a broader filter when com- 
puting the EACF. Again, the analysis is not conclusive for one 
mixed mode at low frequency, but it is able to recover the I = 1 
ridge. We note that the even ridge is affected by the proximity 
of mixed modes. As opposed to the case of KIC 10273246, the 
EACF does not provide any further modes. 

4.8. Rotational splitting and inclination 

Stellar rotation removes the (21 + l)-fold degeneracy of the fre- 
quencies of non-radial modes, allowing for a direct measure- 
ment of the stellar angular velocity averaged over the regions 
probed by these modes, as conveyed by Eq.[4] Using the radii and 
masses computed from model-grid-based methods by Creevey 
et al. (in preparation) together with the estimates of P rot , we 
have computed the ratio of the surface angular velocity to the 
Keplerian break-up velocity, i.e., Q/ s/GM/W, which returned a 
value of approximately 1 % for both stars, indicating that these 
are most likely slow rotators. In view of this and given the pre- 
cision achievable from the spectra, we have thus decided not to 
include any second-order effects on the rotational splitting. 

The overall profile of a non-radial multiplet thus consists 
of the sum of 21 + 1 Lorentzian profiles regularly spaced 



Fig. 13. Similar to Fig. [L2]but for the case of KIC 10920273. 

in frequ ency, and sca l ed in height according to the S'i m (i) 
factors dDziembowskil 119771; iDziembowski & Goodel j 1985; 
iGizon & Solankill2003l) : 

^)-P^[P>osi)}\ (15) 
(I + \m\)\ L ' J 

where i is the inclination angle between the direction of the stel- 
lar rotation axis and the line of sight, and P"'(x) are the associ- 
ated Legendre functions. Note that Yum &im(i) — 1, meaning that 
the $\ m (i) factors represent the relative power contained in the 
modes within a multiplet. 

While we are not able to robustly constrain the rotational 
splitting and inclination for both stars, we are however in a po- 
sition to impose loose constraints on these parameters. Figures 
[14] and [15] map the two-dimensional posterior probability distri- 
butions of these parameters respectively for KIC 10273246 and 
KIC 10920273, based on the samples from a MCMC analysis 
of the ten-month-long time series by IAS_OB. We have overlaid 
each of these correlation maps with curves representing the esti- 
mate of P mt given in Sect. l4.2l and the P ro t(0 relation of Creevey 
et al. (in preparation), obtained by combining the projected rota- 
tional velocity (v sin z) with the stellar radius. 

For KIC 10273246, a comparison of the estimate of the rota- 
tional period with the P ro t(0 relation implies that i > 20°, which is 
also corroborated by the underlying correlation map. Moreover, 
we notice that the marginal posterior probability distribution of 
v ro t is unimodal and points toward an interior rotating roughly as 
fast as the surface. We are however cautious not to claim to have 
robustly constrained the rotational splitting, since changes in the 
model - as well as in the priors on its parameters - used to rep- 
resent the background signal may lead to significant alterations 
in the correlation map. 

For KIC 10920273, no constraints on i are possible from 
a comparison of the estimate of the rotational period with the 
P lot (0 relation. We should nevertheless remind the reader that the 
estimate of P rot is somewhat doubtful (see Sect. I4.21 i. Although 
not completely apparent from Fig. [15] the marginal posterior 
probability distribution of i favors the scenario of a star seen 
pole-on (i.e., i w 0°). Consequently, according to Eq. [15] this 
would make inviable an inference of the rotational splitting. 
Going back to the discussion in Sect. 14.21 the possibility that this 
star is seen pole-on gains strength in explaining the low-SNR 
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Fig. 14. Two-dimensional posterior probability distribution of 
y rot and i for KIC 10273246 (the colorbar decodes the probability 
density level). The prior on i is uniform over the plotted range, 
whereas the prior on v rot is uniform over the interval 0-2 yuHz 
with a decaying Gaussian wing for higher values of the splitting. 
The dot-dashed line marks the estimate of the stellar rotational 
period given in Sect. 14.21 Furthermore, the P mt {i) relation (solid 
curve) of Creevey et al. (in preparation) is shown together with 
the corresponding confidence interval (dotted curves). 




20 40 60 80 

nclination (°) 



Fig. 15. Similar to Fig. QUbut for the case of KIC 10920273. 
Given the possibility that vsin/ * kms -1 , the lower bound of 
the confidence interval of P ro t(0 will coincide with the horizontal 
axis. 

peaks at the low-frequency end of the PDS. There is however an 
alternative interpretation of the correlation map shown in Fig. [15] 
that should not be neglected: the high sensitivity of the parame- 
ter i to rea lization noise may lead to what is known a s fit locking 
at 0° (e.g. JGizon & Solankill2003HBallot et alJl2008h . 

5. Summary and conclusions 

The evolved Sun-like stars KIC 10273246 and KIC 10920273, 
respectively of spectral types F and G, were observed with 
Kepler at short cadence for approximately ten months (from Q0 
to Q4) with a duty cycle in excess of 90%. The light curves used 
in the present analysis have been corrected for instrumental ef- 



from multi-month Kepler photometry 

fects in a manner independent of the Kepler Science Pipeline. 
Both stars are relatively faint and display low SNR in the p-mode 
peaks. 

Different fitting strategies have been employed to extract 
estimates of p-mode frequencies as well as of other individ- 
ual mode parameters, from which we have selected frequency 
lists that will help constraining stellar models. A total of 30 
and 21 modes of degree I = 0,1,2 have been identified for 
KIC 10273246 and KIC 10920273, respectively. These span at 
least eight radial orders. Furthermore, two avoided crossings 
(I- 1 ridge) have been identified for KIC 10273246, whereas one 
avoided crossing plus another likely one have been identified for 
KIC 10920273. These avoided crossings yield strong constraints 
- although model-dependent - on stellar age. Such results con- 
firm previous expectations that asteroseismology of solar-type 
KASC survey targets is possible down to apparent magnitudes 
of 1 1 and fainter, provided w e work with a multi-month time 
series (e.g.,[Stello et a ll2009bh . 

The peak-bagging results presented in Sect. 14.51 are based 
on the analysis of data from the first seven months of observa- 
tions. A smaller number of individual fitters also analysed the 
ten-month-long time series (from Q0 to Q4). In spite of real- 
ization noise being expected to scale as 1 / V r (Libbrecht 1992; 
iToutain & Appourchau"xlll994t IStahn & Gizonll2008l) thus mak- 
ing it possible to increase measurement precisions, there was 
barely any gain in terms of the number of modes detected. 
This is actually not surprising as, once a mode is resolved, 
the background-to-signal rat io in the power sp ectrum cannot 
be improved with time (e.p.. lChaplin et akll2003l) . The only ex- 
ceptions would be (i) a radial mode at 690.34 ± 0.23 pHz for 
KIC 10273246, which is, however, detected with the EACF and 
shown in Fig. [12] and (ii) a I = 2 mode at 8 18. 73 ±0.23 pMz for 
KIC 10920273. 

Good agreement is found between the observed and pre- 
dicted mode amplitudes for the F-t ype star KIC 10273246, base d 
on the revised scaling relation of Kieldsen &Beddind d2011l) . 
This is a particularly interesting result that calls for further tests 
of this scaling relation using the large sample of Kepler stars. 

Despite blending of the multiplet componen ts of non-radial 
modes, i.e., (F) > v rot (e.g., Bal lot et ail 12006). we believe to 
be possible to impose loose constraints on the rotational split- 
ting and stellar inclination for both stars. These constraints are 
based on a combined analysis involving correlation maps, P ro t(0 
curves and estimates of P mt . 

The results presented here point towards KIC 10273246 and 
KIC 10920273 being most likely evolved main-sequence stars. 
The global asteroseismic parameters reported for these stars, to- 
gether with a detailed atmospheric analysis, should allow con- 
straining their radius, mass and age with considerable precision 
(Creevey et al., in preparation). Further insight into the physics 
of these evolved solar-type stars - based on detailed modeling 
and inversion techniques - is now possible due to the high qual- 
ity of the seismic parameters found. 
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Appendix A: Implementing Peirce's criterion 

Peirce's criterion is an exact rule for the rejection of doubt- 
ful observations derived from the fun damental principles of 
Probability theory. Ouoting rPeircd (1 1 8521) : "The proposed obser- 
vations should be rejected when the probability of the system of 
errors obtained by retaining them is less than that of the system 
of errors obtained by their rejection multiplied by the probability 
of making so many, and no more, abnormal observations." 

Logic calls for an iterative assessment of the rejection when 
one or more observations are rejected. The iteration stops when 
no fur ther improvement is possible. Based on the work of lGouldl 
(118551) . we have implemented Peirce's criterion as follows: 



1. 



Compute the mean, x, and the standard deviation, <x, for the 
observational sample {x,}; 

Compute the rejection factor r from iGouldl f 1 8 5 5 ) assuming 

one doubtful observation; 

Reject observations satisfying \x, - x\> ro", 

If n observations are rejected then compute a new rejection 

factor r assuming n + 1 doubtful observations; 

Repeat steps 3 to 4 until the number of rejected observations 

no longer increases. 
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Table 7. Linewidths, heights, and bolometric amplitudes of radial modes for KIC 10273246 returned by the best fit (these are not 
provided for the radial mode of highest frequency since convergence was not properly achieved). 



Frequency 


Linewidth 


Uncertainty 


Height 


Uncertainty 


Amplitude 


Uncertainty 


0"Hz) 


OuHz) 


0"Hz) 


(ppm 2 yuHz~') 


(ppm 2 /iHz~') 


(ppm) 


(ppm) 


737.90 


1.6 


+0.5/-0.4 


10.4 


+3.3/-2.5 


5.48 


±0.81 


785.40 


1.7 


+0.4/-0.3 


11.4 


+3.0/-2.4 


5.90 


±0.80 


833.90 


1.5 


+0.3/-0.3 


19.2 


+4J/-3.8 


7.17 


±0.75 


883.50 


1.8 


+0.4/-0.3 


15.9 


+3.6/-2.9 


7.14 


±0.75 


932.70 


2.3 


+0.5/-0.4 


10.6 


+2.4/-2.0 


6.60 


±0.77 


981.10 


1.9 


+0.5/-0.4 


12.9 


+3.3/-X6 


6.61 


±0.73 


1030.70 


2.6 


+0.8/-0.6 


5.5 


+ 1.5/- 1.2 


5.09 


±0.83 



1079.30 



Table 8. Linewidths, heights, and bolometric amplitudes of radial modes for KIC 10920273 returned by the best fit. 



Frequency 


Linewidth 


Uncertainty 


Height 


Uncertainty 


Amplitude 


Uncertainty 


0"Hz) 


(^Hz) 


0"Hz) 


(ppm 2 yuHz _1 ) 


(ppm 2 yuHz~') 


(ppm) 


(ppm) 


826.66 


0.88 


+0.18/-0.15 


9.10 


+4.27/-2.91 


3.75 


±1.32 


882.77 


0.89 


+0.16/-0.14 


18.56 


+5.30/-4.12 


5.18 


±1.08 


939.58 


0.78 


+0.13/-0.11 


26.28 


+6.91/-5.47 


5.74 


±1.02 


997.14 


0.86 


+0.14/-0.12 


42.21 


+9.71/-7.90 


7.56 


±0.99 


1054.33 


1.05 


+0.20/-0.16 


15.43 


+4.57/-3.53 


5.13 


±1.11 


1111.51 


0.97 


+0.20/-0.16 


12.83 


+4.36/-3.25 


4.54 


±1.16 


1170.77 


0.91 


+0.20/-0.17 


3.77 


+3.26/- 1.75 


2.97 


±1.88 


1226.34 


0.92 


+0.20/-0.17 


5.99 


+3.34/-2.14 


3.26 


±1.46 
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